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I. INTRODUCTION 



Light front Hamiltonian methods have been apphed successfully to theories such as 
Yukawa theory, QED, 0^ theory and (l+l)-dimensional QCD to obtain the mass spec- 
trum and wave functions However, there are still unsolved issues which must be 
focused upon before one can develop non-perturbative methods for light-front quantum 
chromodynamics (LFQCD)[5|. One of these important issues is the problem of infrared (IR) 
divergences 5|, ISj. In quantum electrodynamics, infrared divergences are eliminated if one 
uses appropriate initial and final states of charged particles with a suitable superposition of 
an infinite number of photons. The issue of cancellation of IR divergences at amplitude level 
has been addressed by Chung j3| and Kulish and Faddeev (KF) . KF proposed a method 
of asymptotic dynamics wherein one replaces the free Hamiltonian by an asymptotic Hamil- 
tonian that takes into account the long range interaction between incoming and outgoing 
states. This asymptotic Hamiltonian is used to construct an asymptotic evolution operator 



T exp 



I / Vas{t)dt 
' TOO 



(1) 



and thereby a set of coherent states 



\n; ±) = Qj^\n) 



(2) 



It was shown that the transition matrix elements formed using these asymptotic states are 
free of IR divergences. Subsequently, KF method was discussed in the context of QCD by 



various authors |9Nl3|. 

A coherent state formalism in light-front field theory (LFFT) has been discussed by 

nn 

various authors [1^, 115|, in particular, in the context of cancellation of IR divergences in 



QED and QCD at lowest order 



IGMISI]. One of us proposed the use of coherent state 



formalism in LFFT due to its usefulness in establishing the cancellation/ non-cancellation 



of IR divergences in Hamiltonian formalism 16 



A coherent state formalism for LFQED 



in LF gauge was developed and applied to show the cancellation of IR divergence in one loop 
vertex correction 1^ . Recently, we have extended the formalism to demonstrate cancellation 
of IR divergences up to O(e^) in fermion self energy in Ref. 2^, henceforth referred to as X. 



LF quantization is performed usually in light 



tages when applied to non-abelian theory 



21 



ront gauge A'^ = due to its many advan- 
22| , in particular, the absence of ghost fields. 
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One of the reasons for using LF gauge is that one is able to solve the constraint equation for 
the non-dynamical part of the fermion field. Gauge independence of LFQED calculations 
has been addressed in literature. In particular, gauge dependence of non perturbative cal- 
culations has been studied and gauge invariance of physical quantities has been confirmed 



23|. In Ref. 



24j . mass eigenvalue problem in QED is discussed in an arbitrary covariant 



gauge. In this work, the theory is regulated by introducing PV photons and electrons. It has 
been shown that in Feynman gauge, one PV electron and one PV photon are sufficient for 
canceling the instantaneous fermion interactions and the A"*" dependent terms cancel from 



the constraint equation 2J] . 

In X, we have shown the cancellation of IR divergences up to 0(6^) in LF gauge using 
coherent state basis. Here, we extend our analysis to show this cancellation in Feynman 
gauge. We calculate fermion self energy in LFQED up to O(e^) in Feynman gauge. We 
show that in Feynman gauge also the true IR divergences in lepton mass renormalization 
are cancelled up to O(e^) if one uses coherent state basis for evaluating the transition matrix 
elements. The plan of the paper is as follows: In Section II, we present the Hamiltonian 
of PV regularized LFQED in general covariant gauge 24l| and calculate the O(e^) fermion 
mass renormalization using light-cone-time-ordered perturbation theory(LCTOPT) in Fock 
basis. We demonstrate the appearance of true IR divergences in the form of vanishing light- 
cone energy denominators. In Section III, we obtain the form of coherent states from this 
Hamiltonian using the method of asymptotic dynamics. In Section IV, we calculate dm^ 
in lowest order using the coherent state basis and show that the extra contributions due 
to emission and absorption of soft photons indeed cancel the IR divergences in dm^. In 
Section V, we calculate dm? up to 0{e^) in Fock basis and identify the true IR divergences 
in it. In Section VI, we perform the same calculation in coherent state basis and show the 
cancellation of IR divergences. Section VII contains a summary and discussion of our results. 
Appendix A and B contain the details of the calculation of transition matrix elements in 
Fock basis and coherent state basis respectively. 
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II. PRELIMINARIES 



A. Light-Front QED in Feynman gauge 



We begin with the QED Lagrangian in Lorentz gauge with an arbitrary gauge parameter 
( and additional PV fields 



2|: 



1=0 

2 



(3) 



i=0 



Here 



(4) 



1=0 



i=0 



i = corresponds to the physical field and i = 1 and 2, to PV fields. The fermion and 
photon fields have masses mj and /Xj respectively and /Uq = for the physical photon. /3j and 
are coupling coefficients which satisfy 

2 2 

Y.^-m = o, 5^(-im = o. (5) 

i=0 i=0 

The light-front Hamiltonian density, for a free massive vector field with mass is given by 



n = n\c=i + ^(1 -0{d-A){d-A- 2d^A+ + 2d^A^), 
where the first term of the Eq. ([H]) gives the Feynman gauge Hamiltonian density 



(6) 



(7) 



At=0 



Thus, the light-front QED Hamiltonian in Feynman gauge = 1) is 



P- = H = Ho + V , 



(8) 



where the free Hamiltonian is given by 



1=0 s 
2 3 

Z=0 A=0 



2p+ ' 



s) + s)cii(p, s)) 



(9) 
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and V has the form 



V = e f (f^^dx- I [dp][(ip\[dk] J2 Yl s') + e-'P-'^Wiip, s')d,{p, s')] 

s,s',\ 

X7^[e-^^-^^,(p, s)b,ip, s) + e'^-'v,{p, s)d]{p, s)]el{k)[e-^''-^ai{k, A) + e'^-^a\{k, A)], (10) 



where 



[dp] 



dp'^ 



(11) 



> (2^)1 Jo ^ 

Summation over i, j, / runs from to 2. Since PV- field contributions have been added to the 
Lagrangian, the instantaneous fermion terms cancel 2S\- The instantaneous photon terms 
that appear in light-cone gauge are also not present in Feynman gauge. The polarization 
vectors e^'*''' have an additional flavor index because they depend on the mass of the photon 
flavor. 



B. Lepton Mass Renormalization in Light-Front QED 



In light-front time ordered perturbation theory, the transition matrix is given by the 
perturbative expansion 



T = V + V- 



1 



-V + - 



[12) 



p--Ho 

The lepton mass shift is obtained by calculating Tpp which is the matrix element of the 
above series between the initial and the final lepton states \p, s) and is given by 



We expand Tpp in powers of as 



Sm^ = p+ ^ Tpp 



Tpp = T^^) + r(^) + ■ ■ • 



(13) 



(14) 



In general, T*^"^ gives the 0(e^") contribution to lepton self energy correction. Here, the 
initial (or final) lepton momentum is 



p 



' 2p-^ 



Momentum of internal photon line 

k 



/ = 0,1,2 



(15) 



(16) 
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Momentum of the internal fermion line is 



Pi 



P , 



2p+ '' 



P± 



i = 0,1,2 



O(e^) correction is obtained from 



T«=T«(p,p) = (p,s|\/ 



■V\p,s) 



(17) 



(18) 



Note that 



T«=TW(p,p)=Ti, + Ti 



lb 



(19) 



where Tia and Tn, are O(e^) contributions from standard three point vertices involving 
physical photon and PV photons respectively and are represented by the diagrams in Fig. 
1(a) and Fig. 1(b). The wavy arc represents physical photon while the spiral arc represents 




(a) 




FIG. 1: Diagrams for O(e^) self energy correction in fock basis corresponding to Tia and Tif,. In 
Fig. (a), wavy line corresponds to physical photon (i = 0) and in Fig. (b) curly line corresponds 
to PV photon (i = 1,2). 



both PV photons {j = 1,2). The index i on the internal fermion line takes values i = 0, 1 
or 2, where i = corresponds to the physical fermion field while i = 1 and 2 represent PV 
fermion fields. We are interested in the true IR divergences which arise due to vanishing 



energy denominators in TOPT 



16|. 



The energy denominator in Tu, is 



p — {p — k 



PI 



'i + {p± - h^f + + ^2 ^ ^2 



1 i 



2p^ 



2{p+ - kt) 



2k+ 



p ■ ki — + mf 



p+ - kt 



(20) 



One can show that this energy difference is zero only when rrii = m. Therefore, up to 
0(6^) the energy denominator cannot be divergent if flavor changing vertices are involved. 
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Moreover, due to non-zero mass of photon PV field, in the limit kf — )■ 0, ki^ 0, p- ki -/^ 
irrespective of the flavor of fermion field. Thus Tif, cannot have such IR divergences and we 
do not need this diagram for our discussion. Neglecting Ti^, O(e^) transition matrix element 
contributing to fermion self energy reduces to 

Tia{p,p) = {p,s\Vi ^ Vi\p,s) (21) 

P - Hq 

Here Vi is same as in X i.e. 3-point vertex involving physical fermion and physical photon 
(i = 0) only. Proceeding as in X, 



2(2jr)3i '"i k+p* Mp--p--K 
where pi = p — ki. 

In the asymptotic limit, kf — )■ 0, ki^ — > 0, we obtain 

It is to be noted that the diagrams with PV photon fields and/ or flavor changing vertices do 
not contribute to IR divergences, as in both the cases the energy denominator is necessarily 
non-zero. Hence, we need to consider i = case only for Fig 1(a). The denominator vanishes 



in the limit k^ — )■ 0, kix — ?■ i .e p ■ fci — )■ 0, thus leading to IR divergences 16|, 120 1 



It was pointed out in Ref. [2^ that in the infinite PV mass limit a tree level diagram 
involving an intermediate PV fermion of mass mi reduces to four point instantaneous vertex 
as illustrated in Fig. 2. Similarly, Fig. 1(a) with i = 1, reduces to the instantaneous diagram 




mi — > OD 



FIG. 2: In the infinite PV mass limit the PV fermion line reduces to an instantaneous four point 
interaction term denoted by a dash on the fermion line . 

in Fig. 1(b) of X in the infinite PV fermion mass limit. However, such diagrams do not 
contain a vanishing energy denominator and hence are not relevant for the present discussion. 
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III. COHERENT STATE FORMALISM AND INFRARED DIVERGENCES 



In X, it was shown that the true IR divergences in self energy correction up to O(e^) get 
cancelled if one uses coherent state basis in LFQED in LF gauge. We will prove the same 
result in Feynman gauge here. For that purpose, we will now obtain the form of coherent 
states for LFQED Jn Feynman gauge by the method used in Ref. js] for equal time theory. 



As shown in Ref. 2J], when one writes the Hamiltonian in terms of independent degrees 
of freedom, the non-local terms do not appear and hence, in this formalism, there are no 
instantaneous interaction terms in LF Hamiltonian. 



The light-front time dependence of the interaction Hamiltonian is given by 16 1 



4 r 

V{x^) = eJ2 diy^[e-""^^^K{iya) + e'"''^^ hUiy^)] (24) 

a=l ^ 



hai^a) are three point interaction vertices. For example, 

2 

^1 = E E ^)«'(^' ^)^(^' ^')7^«,b, , (25) 

i,j,l=0 sX 

and Ua is the light-front energy transferred at the vertex ha- For example, 

-i = p7 + ^r-pr = ^|^ (26) 

is the energy transfer at eej vertex. The integration measure is given by 

1 

(2^) 

and being fixed at each vertex by momentum conservation. 

In the limits \x~^\ — )■ oo, non-zero contributions to V{x^) come only from regions where 
Pa 0. It is easy to see that U2 and z/3 are always non-zero and therefore, h2 and 
do not appear in the asymptotic Hamiltonian. As shown earlier in Eq.( l20|) the energy 
differences ui and cannot be zero at flavor changing vertices or at vertices involving 
massive photon fields. This means and z/4 can be zero only for the physical photons 
and fermions but not for PV photons and fermions. Thus, the IR divergent contribution to 
asymptotic Hamiltonian comes only from the terms with i = j = I = in Eq. (Elj). Thus, 



d^=^r^l^^ (27) 



the 3-point asymptotic Hamiltonian is defined by the following expression 16 1 



VUx+) = eY, [ c?^«0A(A:)[e-*^"^+/if (^«) + e'"''''^ U^^^u^)] (28) 



where 0a (fc) is a function which takes value 1 in the asymptotic region and is zero elsewhere. 



The superscript on ha\va) indicates the z = j = / = part of V{x^). Eq. fl28|) is same 



as V„Jx^) in Ref. 



161 and 3-point vertex part o: 



for obtaining Vas{x^) can be found in the Refs. 
from the asymptotic Hamiltonian: 



16|, 



Vas{x^) in X. The detailed procedure 



20|. The asymptotic states are obtained 



|n; ±) = VL^\n) 



(29) 



where is the asymptotic evolution operator and \n) is the Fock state, fi^ is defined by 



= T exp 







« / Vas{x^)dx^ 
I =Foo 



(30) 



Following the procedure in X, the asymptotic states are found to be 



VL^\n) =exp 



-e J dpU'p^ J 



- (27r)3/2 7 ^ 



dk^ 



Here 



[/(fc, A : p)a\k, A) - r (fc, A : p)a{k, A)]p(p) 



p{p) = J2ibKpMp)-dl{p)d,{p)), 

i=0 



\n) 



and 



fik,X:p) 



a{k, A) = ao(fc, A) 
p,6^(fc) /fc+A 



-k ^ 



(31) 

(32) 
(33) 

(34) 



p ■ k \ p~^ y V 

One must notice that the form of asymptotic state is simpler here as compared to that in 
LF gauge (discussed in X) as 4-point instantaneous interaction does not appear in Feynman 
gauge Hamiltonian. Furthermore, there is no contribution to asymptotic Hamiltonian from 
PV fields as they are massive and hence p~ — k~ — {p — k)' in this case. One can use 
Eqs. 02]), ( pTj) and ( 15^ to obtain the form of coherent states 

d^ki r dk^ 



=exp 



-.[fik,X,p)a\k,X) - r{k,X,p)a{k,X)] 



\n) (35) 
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\ fcl fel / 



\ 

^ 



(P,s) (p,s) 

FIG. 3: Additional diagrams in coherent state basis for O(e^) self energy correction corresponding 
to T' 

IV. LEPTON MASS RENORMALIZATION UP TO ©(e^) IN COHERENT STATE 
BASIS 

The self energy contribution up to O(e^) in coherent state basis is given 
where T^^^ is defined in Eg. (1181) and T'^^^ arises from O(e^) term in 

T'ip,p)={p,s: f{p)\V\p,s: f{p)) (36) 

and is represented by Fig. 3. A soft photon in coherent state is shown by the dotted line in 
the Feynman diagram in Fig. 3. The detailed calculation has been discussed in X. Here we 
present only the result: 



where the prime indicates the correction due to additional terms in coherent state basis. 
The energy denominator in Eq. ( 1371) vanishes in the limit ^ 0, ki± — thus leading to 
IR divergences. Adding Eqs. ( l23l) and ( 137|) . these true IR divergences get cancelled and the 
O(e^) lepton mass correction is IR divergence free. 

V. LEPTON MASS RENORMALIZATION UPTO 0{e'^) IN FOCK BASIS 

We will now calculate O(e^) lepton mass correction in Fock basis. Transition matrix 
element for O(e^) correction to self energy is given by 

T(2) = (38) 

where 

p - Hq p - p - 

(40) 



2 



dkt{p-e{h)yQ^{h) 
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Note that due to presence of PV fields, we get an additional set of diagrams involving 3- 
point vertices in addition to diagrams in X. All these diagrams can be evaluated in the 
standard manner by inserting appropriate number of complete sets of intermediate states. 
Moreover, in present formalism, the four point instantaneous interaction term are absent in 
the Hamiltonian and therefore diagrams in Figs. 4, 5 and 6 in X are not present here. We 
will not reevaluate the diagrams already evaluated in X and discuss here only the additional 
diagrams which appear due to introduction of PV fields. We give the details of the calcu- 
lation in Appendix A and present here only the results for the new diagrams. Additional 




k {p.. s) 



FIG. 4: Additional diagrams for O(e^) self energy correction in fock basis corresponding to 
(Complete set of diagrams consists of these and Figure 3 in I) 

contribution to due to PV fields is given by 



rr ^ Trip, P) = + + + 



Ad 



(41) 



where , , and correspond to Figs. 4(a)-(d) and are given by Eqs. (lA3p . 
(IA6p and f lASp . Using expressions for energy denominator from Eqs. flAip and flA2p we 
obtain 



{6 m 



2\PV 
Ua 



2(2^)6 7 kt ktpfptpt 



32{p -Pk-k^){p -Pi -k^){p -p. -k^ 



(42) 
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2(27r)6 7 kt Hpfpjpt 

32 (p- -Pk - ^2~)(P~ - P7 - ^r)(P" - Pj - K - k^) 

(43) 



kj k+pjpjpl 

Tr[4^^{k2){A + mk) 4^^{k2){^'j + mj) (A;i)( + m,) {k^) { ^ + m)] 
32{p~ -pl - k2){p-^ - Pi - kl){p'- - p'r) 

(44) 

where p'j = p and pi, pj and pk have been defined in Eqs. (]A4p . (lASP and (]A7p respectively. 
The last term is IR convergent, as the energy denominator involved can never be 
zero, and hence it is not needed for our discussion. Note that the divergence structure of 
these diagrams and those in X is different i.e ((5m^)£j^, (5m^)ff,^ and {6m^)f^ can have IR 
divergences only when p ■ ki ^ i.e kf ^ 0, kij_ -^0, i = j = k = and / = i.e. only 
when the diagram involves physical fermions and photons. Also, as the PV field is massive, 
the energy denominator involving its k~ momentum cannot be zero and therefore we need 
not consider the limit A;^ — 0, k2_L — ?■ and the combined limit kf ^ 0, ki_L — )■ 0, A;^ — )■ 
0, k2_L — as we did in X. As shown in Appendix A, the IR contribution from diagrams in 
Figs. 4(a) and 4(b) is given by 



e 



4b 



rf^ki^d^k^x / dkt [2{p ■ e{k,))\p ■ ejk^))' - {p ■ k2){p ■ e{k,)y, ^^^^ 



{2nf J J kt kt A{p ■ ki)\p ■ k2) 

and the contribution from Fig. 4(c) is given by 

2^Py /,2, ,2, f dkt dkt 2{p . e{k,))\p . e{k2)f - {p ■ k2){p ■ e{k,)Y 



(27r)6 J J kt kt 8p+ 

(46) 



P ^ Pk 



{p-kiY{p-k2) {p ■ ki){p ■ k2Y _ 
Here we have used Heitler method {29 1 discussed in Appendix A to deal with the vanishing 
denominator {p~ —p'~). 

In addition to diagrams in Fig. 4, there are O(e^) diagrams containing only physical 
photons but PV fermions in the intermediate states. These will appear the same as diagrams 

12 



in Fig. 3(a) of X except that the internal fermion hnes will correspond to PV fermions. One 
can draw these diagrams by replacing the curly lines by wavy lines in Fig. 4(a)- (c) here. 
Such diagrams will not involve vanishing energy denominators due to the presence of flavor 
changing interactions as discussed in Section 11. However, in the infinite PV fermion mass 
limit, such diagrams will reduce to diagrams involving instantaneous interaction vertex. 
This is illustrated schematically in Fig. 5, where Fig. 3(a) of X with the intermediate term 
involving PV fermion field reduces to Fig. 4(a) of X in infinite PV fermion mass limit. 




j k 

3(a) 4(a) 



FIG. 5: In the infinite PV mass limit the diagram on the left reduces to a diagram involving 
instantaneous interaction. Here i = 1 while j = k = 0. This figure involves intermediate PV 
fermion also while Fig. 3(a) of X involved only the physical fermion terms. 



VI. LEPTON MASS RENORMALIZATION IN COHERENT STATE BASIS UP 

TO 0{e^) 

In this section, we calculate the 0{e^) lepton mass correction using coherent state basis. 
We will show that the IR divergent contribution from the additional diagrams in coherent 
state basis exactly cancel the IR divergences arising due to vanishing energy denominators 
calculated in Section V. In coherent state basis, 0{e'^) correction to self energy is given by 

y(2) ^ j.^ 

where is O(e^) term in {p, s: f{p)\V ^_^^^ V ^_^^^ V\p, s: f{p)) represented by Fig. 6 here 
and Fig. 4 in X. We present the details of calculation in Appendix B and give below only 
the contribution of Fig. 6 to [drn?)' . Cancellation of divergences in diagrams involving only 
massless photon has been given in X and we will not repeat it here. The contribution of Fig. 
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(a) 




{P,s) {p,s) 



FIG. 6: Additional diagrams in coherent state basis for O(e^) self energy correction corresponding 
to 



6 is given by 

= {5m%^ + {5m%, + {5m%^ (47) 

((5m^)g^, (5m^)g^ and (^m^)g^ have been evaluated in Appendix B and we give below only 
the result, 

f J2i, j2i, f ^^i" '^^2' 



-QAih) (48) 



1 2 

[2(p ■ 6(A;i))^(p ■ e{h)Y - (p ■ fc2)(p • e{h)f 



(W)gft=-^^ J (fki^(fk2± J 



2 

21 



-(*^»-' B4h) (49) 



4(p ■ fci)(p ■ k2)[{p ■ h) + {p ■ h) - ih ■ k2)] 
Adding Eqs. (jH]) and we obtain 

(W)L + (W);,=^ / d^k,,rf^k,, I 

[2(p-6(A;i))2(p-e(fc2))^-b-^2)(p-6(A:i))^] 



A{p-hy{p-k2) 



QAih) (50) 



14 



Also 



(27r)6 



+ 



dkt dkt [2{p ■ e{h))\p ■ e{k2)f - (p ■ k^){p ■ e{k,)f 



{p-kiY{p-k2) {,p-k2Y{p-ki) 



QA{ki 



(51) 



Adding Eqs. ^ and dSO]) we see that {5m^)'[^ + (5m2)f^^ + {5m%^ + (5m2)[j^ is IR finite. 



Similarly, adding Eqs. f l46p and (151 p we find {5m?)^Y + (^"^^)6c ■'■-^ finite. Adding all the 
contributions coming from Figs. 4 and 6 and combining with the results in X, we find that 
the self energy correction up to O(e^) is IR finite. This completes the proof of cancellation 
of true IR divergences up to 0{e^) for fermion self energy correction in coherent state basis 
in Feynman gauge also. 



VII. CONCLUSION 



In X, we have calculated lepton self energy correction in light-front QED in LC gauge up 
to O(e^) and have shown that the true IR divergences get cancelled when coherent state basis 
is used to calculate the matrix elements. In this work, we have obtained the same result 
in Feynman gauge. One can notice that the proof of cancellation is simpler in Feynman 
gauge. The proof can be generalized to general covariant gauges. We plan to address this 
in a future work. 

The cancellation of IR divergences between real and virtual processes is known to hold in 
equal-time QED to all orders. This cancellation was also shown by Kulish and Faddeev 8| 
using the coherent state formalism. KF method leads to cancellation of IR divergences in 



QED to all orders. However, it is well known that the Bloch-Nordseick theorem [25[ does not 
hold in QCD and therefore, in QCD one does not expect to construct an all order proof of 
cancellation of IR divergences along the lines of KF method. Basically, the non-cancellation 
of IR divergences in QCD arises due to the fact that asymptotic states here are bound 
states of quarks and anti-quarks and therefore the asymptotic Hamiltonian to be used in 
KF method should contain the confining potential and should not be just the asymptotic 
Hamiltonian of QCD. An "improved" method of asymptotic dynamics has been introduced 



by McMuUan etal 



26l428l| which takes into account the separation of particles also. The 



improved method has also been discussed in the context of LFQED and LFQCD[1 
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An improved coherent state method in LFFT may be useful from the point of view of 
extracting information about the artificial confining potential which is needed in LF bound 
state calculations jsl- If we use appropriate Hamiltonian of bound states as the asymptotic 
Hamiltonian and develop a coherent state approach based on it, then this approach would 
lead to cancellation of IR divergences in QCD as well. It will be interesting to understand 
this connection between cancellation/non-cancellation of IR divergences and the form of 
asymptotic Hamiltonian. Thus we may be able to get some perception of the form of 
artificial confining potential mentioned by Wilson et. al. in Ref. jsl by understanding the 
structure of IR divergences in LFQCD. 
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Appendix A: Transition matrix element for self energy in fock basis 
We will need the following expressions for energy denominators: 

p--K-{p-k,r = -^^ (Al) 

~ i~ ( 1 p- ki+p- k2- ki- k2 

P -k^ -k^ - [p-ki- k2) = -+ -+ (A2) 

p+ — k^ — ^2 



We now calculate T4 which is defined by Eq. fl39|) and corresponds to Fig. 4. Inserting 
complete sets of intermediate states in T^^, we obtain 

where 

rppv _ / c?^ki^(i2k2_L f dk^dk^ 

^'^ "(2^7 2^ J 32ktktptplpt 

u{p,s)[4^'{ki)i^k + mk) 4^^{k2){^j + mj) 4^'{k2){jf)i + mi) /^'{ki)]u{p, s) ^^^^ 
{P" -Pk - K){p~ - Pi - K)iP~ - Pj - K - ^2") 

with 

Pi=p-ki, (A4) 
Pj = p - ki - k2. (A5) 
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Note that, here pi = pk 
Similarly, 

py f (i^ki^(i^k2_L f dk^dk2 

■J- / 



(2^)6 y 2p+ J 32ktktptpjpt 

(A6) 



u{p,s)[4^^{k2){^k + mk) 4^'{ki){^j+m,) 4^'{k2){h + rui) 4^^{ki)]u{p, s) 



{P~ -Pk - h -Pi -K )ip~ - Pj -k^ -k^) 
with 

Pk=p-k2 (AT) 
py f cPki±cPk2± f dkidk2 



T — 

(2^)6 y 2p+ J 32ktktptptpt 
u{p,s)[4^'{k2){h + mk) 4^'{k2){^'j + m) /^^(A:i)(i>, + m,) 4^^{k,)]u{p, s) 

{p~ - Pt - K){p~ - p'f){p~ - Pk - ^2) 

In the limit, /c]^ — 0, ki^ and i = j = A; = 0, we get the IR divergent contribution. 
In this limit, Eqs. flASP and flA6P can be added such that the denominator reduces to 
{p-kif{p-k2) 

2\pv , /r_2\py f j2i, j2i, f dktdk2 



{Sm')Z + =^ / d^k,^d^k2^ I 



Tr[/^^{ki){^ + m) /^''{k2){h + m) /^^{k2){^ + m) /^'{ki){^ + m)] 



{p ■ kiY{p ■ k2) 



(A9) 



After calculating the trace, Eq. ( ]A9P leads to 

- j2^^ J J It ^ <p-hnp-k2) 

Similarly, leads to 



2{2nfJ -J 32ktktptptpt 

Tr[4^'{k2){h + m) 4^^ik2)i^'j + m) 4^'{ki){^i + m) 4^'{ki){^ + m)] 



(All) 



{P -P^-K){P -P2 )(P -Pk-h) 

where pi and pk are defined by Eqs. (lA4p and (1A7P and p^- = p. Note that this diagram is one- 
particle reducible, and therefore the energy denominator associated with the single-particle 
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state vanishes. We have used the Heitler method [29| for evaluating all such integrals. Using 



this method, we write 
D 



29 



3Q| 



{P~ - P'j ){P~ -Pi - ^1 ){P~ -Pk - ^2 



j dp' 5{p' —p 



V 



{p' - p'j ){p' -Pt -h){p' ~Pk -hY 

Using the relation between distributions 

-_ -p ) = --S\p' -p ). 

[P -Pj ) 2 

and integrating by parts we obtain 

D=y dp'-s(p'- -p-)^ 



2(p - Pi - h f{P - Pk - h ) 2(p - -K){p -Pk-h f 
Thus, in limit I, Eq. (lAlip becomes 



(A12) 



(A13) 



(A14) 



dkt dk+ 2{p ■ e{h))\p ■ e[k2)f - [p ■ k2){p ■ e{k,)f 



+ 



Pk 



8p^ 



_{p ■ kiY{p ■ k2) ip ■ ki){p ■ k2Y _ 
The traces are calculated using Mathematica. 



(A15) 



Appendix B: Transition matrix element for self energy in coherent state basis 

The contribution corresponding to diagrams in Fig. 6 is given by 



2(2^)6 J "^^ J IGk^k^ 

Tr[4^^{k2){h + m) 4^'iki){^i + m) 4^'{ki){^ + m)]{p ■ e^'{ki))QA{ki) 



ip-kin{p-k,) + {p-k2)-{ki-k2)] 



(Bl) 



{6m 



2\' 



(i^ki_L(i^k2^ 



'^^ 2(27r)6 J " "'^ J 16k+k+ 

Tr[jt^'{k2){h + m) 4^'ih){h + m) 4^'{k2){^ + m)]{p ■ e^^(A;i))9A(fci) 
{p ■ ki){p ■ k2)[{p ■ ki) + {p ■ k2) - {ki ■ k2)\ 



(B2) 
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Calculating the trace, adding Eqs. (IBip and (lB2p and taking limit I, we obtain Eq. ( !50|) . 



{p^ - Pi -h - P'j ){P~ -Pk -^2 



(B3) 



Now using Heitler method 29|, we obtain Eq. ( 15T|) . 
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